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A discretization of a continuum theory with constraints or conserved quantities is called mimetic 

if it mirrors the conserved laws or constraints of the continuum theory at the discrete level. Such 

discretizations have been found useful in continuum mechanics and in electromagnetism. We have 

recently introduced a new technique for discretizing constrained theories. The technique yields 

discretizations that are consistent, in the sense that the constraints and evolution equations can 

be solved simultaneously, but it cannot be considered mimetic since it achieves consistency by 

' determining the Lagrange multipliers. In this paper we would like to show that when applied to 

general relativity linearized around a Minkowski background the technique yields a discretization 

^ , that is mimetic in the traditional sense of the word. We show this using the traditional metric 

, O ■ variables and also the Ashtekar new variables, but in the latter case we restrict ourselves to the 

' Euclidean case. We also argue that there appear to exist conceptual difficulties to the construction 

of a mimetic formulation of the full Einstein equations, and suggest that the new discretization 

scheme can provide an alternative that is nevertheless close in spirit to the traditional mimetic 

■ formulations. 
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' Continuum theories, either mechanical systems or field theories, usually have conservation laws and sometimes 
, constraints. When one discretizes the equations of these theories, for instance in order to solve them numerically on 
a computer, or for "quantization on the lattice" purposes, the resulting discrete equations will usually fail to preserve 
the conserved quantities of the continuum theory upon evolution. Similar comments apply to constraints. Although 
one may have discrete equations resulting from discretizing the constraints of the continuum theory, if one chooses 
initial data that solves these equations exactly, they will fail to be solved upon discrete evolution. 
?-H Mimetic discretizations are discretizations of continuum theories that preserve conserved quantities or constraints 

in the discrete theory that mimic those of the continuum theory. There is quite a body of literature on mimetic 
^ discretizations in the context of continuum mechanics and electromagnetism. The literature on Hamiltonian lattice 
QCD implicitly considers a mimetic discretization of Yang-Mills theory, although this fact is not usually emphasized. 
rS Some authors have considered the question of whether mimetic discretizations of general relativity can be con- 

■ structed i^]. It is well known that if one discretizes the Einstein equations, the Hamiltonian and momentum 
constraint, which should hold for all time if satisfied initially (ignoring for the moment the issue of spatial bound- 
aries), fail to do so in the discrete theory. Although there has been success in generating mimetic formulations of 
linearized relativity, it appears unlikely that something similar will be available for the full theory (or even for the 
linearized theory on non-trivial backgrounds or slicings). This is due to the fact that discretized derivatives fail to 
satisfy Leibnitz' rule and therefore the nonlinear terms when discretized do not have properties that mirror those of 
the continuum Q. 

We have recently introduced a new approach to the discretization of theories, particularly of theories with constraints 

0. 13 called "consistent discretization". The technique guarantees that the resulting discrete equations are compatible, 

1. e., they admit a common set of solutions (something that is not generically true if one discretizes the equations of 
a constrained theory). The technique has been tried out in the context of cosmological solutions of the Einstein 
equations jfl], of BE theory and of Maxwell and Yang-Mills theories on the lattice Q. Current investigations are 
testing it for the Gowdy models. 

In this paper we would like to show that the technique we proposed, when applied to the Einstein theory linearized 
around a Minkowski background yields a discrete formulation that is mimetic. That is, the discretized constraints are 
exactly preserved under evolution without determining the Lagrange multipliers. We first consider linearized general 
relativity in terms of the traditional metric variables. We then consider it in terms of Ashtekar's variables, which 
have the advantage of being closer to the discretizations used in Yang-Mills theories (although in this case we restrict 
ourselves to Euclidean general relativity) . 



2 



In the consistent discretization scheme, equations are discretizcd with variables evaluated at two (or more) different 
levels in time. This includes the constraints of general relativity, which in the discrete theory therefore can only be 
viewed as "pseudo" constraints (we will reserve the word constraint for expressions that involve all variables evaluated 
at the same instant of time, as in traditional canonical terminology). These equations, with variables discretized at 
mixed instants of time are the equations that are solved by the consistent discretization scheme. Of course, if one 
is in a regime in which the time-step is small, then satisfying the pseudo-constraints implies that the usual discrete 
constraints (with all the variables at the same time-step) are approximately satisfied as well. Therefore the resulting 
scheme cannot be strictly called mimetic, although it approximately is. We will show that if one uses a discretization 
for general relativity that is mimetic in the linearized case, one further improves the accuracy with which the consistent 
scheme for the full non-linear theory satisfies the constraints. This encourages further studies of these discretization 
schemes in the context of numerical applications. 

In the next section wc will present a brief summary of the consistent discretizations scheme. In the following two 
sections we apply it to linearized gravity, first with the traditional variables and then with the Ashtekar variables. 
We end with a discussion and proposals for further research. 



II. CONSISTENT DISCRETIZATION OF CONSTRAINED THEORIES 

We illustrate the technique with a mechanical system for simplicity, but there is no problem working it out for field 
theories, since upon discretization the latter become mechanical systems. We assume we start from an action in the 
continuum, written in first-order form. 



J L{q,p)dt 



(1) 
with 

L{q,p)=pq-H{q,p)-X(t>{q,p) (2) 

where A is a Lagrange multiplier and the theory has a single (it is immediate to incorporate several) constraint 
^{q^p) = 0. The discretization of the action yields S = L{n,n+ 1), where 

L{n, n+l)= Pn{qn+i - qn) - eH{qn,Pn) - K(t>{qn,Pn), (3) 

where e — t„+i — tn and we have absorbed an e in the definition of the Lagrange multipliers. 

We will now view the Lagrangian as the generator of a type 1 canonical transformation between the instant n and 
the instant n+1. In ordinary classical mechanics parlance, given a canonical transformation between a canonical pair 
q,p and a new canonical pair Q, P, the generating function of a type 1 canonical transformation is a hmction of (j, Q, 
F{q,Q) and the canonically conjugate momenta are defined by P = dF/dQ, p = dF/dq. In our case we will view 
qn,Pn,Xn and A„+i as "configuration variables" and will assign to each of them a canonically conjugate 

momentum through the canonical transformation. 





dL{n, n + 1) 


(4) 


dqn+i 




pP _ 


dL{n, n+1) 

dpn+l 


(5) 


p\ _ 


dL{n, n + 1) 

dXj^n+l) 


(6) 


pi = 


dL{n, n+1) 
dqn 


(7) 


pp — 

n 


dL{n, n+1) 
dpn 


(8) 


P^ = 

n 


dL{n, n+1) 

9\n) 


(9) 



If one explicitly computes the partial derivatives with the Lagrangian given, one can eliminate the p, and P^ to 
yield a more familiar-looking set of equations, 

p, p, _ dHiq^,P^^,) 90^(9„,P„VJ 

" ~ ' din ""^ din ' 
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qn+i 



d^^iqn,P^^,) 



(10) 



These indeed look like a discrete version of equations for a system with constraints. However, there are important 
differences. First of aU, notice that as an evolution system the equations are implicit. Secondly, if one solves the 
first two equations one obtains and q as functions of the initial data and the Lagrange multipliers. The last 
equation however, will generically not hold. One will have to choose specific values for the Lagrange multipliers at 
each time-step (and if one is dealing with a field theory at each point in space) for all the equations to be solved. 

Notice that there can be particular cases in which the system does not determine the Lagrange multipliers. For 
instance, consider a totally constrained system like general relativity. There the Hamiltonian vanishes. Suppose now 
that the constraint in (|10|l is only a function of qn- Then the evolution equation for qn+i implies that qn+i = qn and 
the constraint is automatically preserved. Therefore the resulting formulation is mimetic in the traditional sense of 
the word, a constraint that is just the discrete version of the continuum constraint is preserved under evolution by 
the discrete evolution equations. A similar situation develops if the constraint is only a function of P*. 

If the Hamiltonian is non-vanishing, and the constraint depends only on P^j^i then the latter is not automatically 
preserved upon evolution, but it cannot be satisfied by choosing the Lagrange multipliers either since they drop out 
from the relevant evolution equations. On the other hand if the constraint is only a function of qn , its preservation 
could be enforced by choosing the Lagrange multipliers (the asymmetry between P'' and q in this treatment comes 
from the fact that we chose to write the equations as "propagating forward" in time, if one had chosen to propagate 
backwards, the roles of q and P^ in this discussion would be reversed). 

Summarizing, the consistent discretization technique consists of discretizing the action and working out the resulting 
equations of motion for the discrete theory from it through the canonical transformation that implements time 
evolution. The resulting evolution equations (and constraints) are made a consistent set of nonlinear algebraic 
equations by considering the Lagrange multipliers as dynamical variables one has to solve for. In particular situations, 
the Lagrange multipliers are not determined by the equations. In such cases the resulting set of equations and 
constraints has to be consistent since it has been derived from a variational principle and the resulting discrete theory 
is mimetic in the traditional sense of the word: the constraints are automatically preserved upon evolution. In the 
other case, when the Lagrange multipliers are determined the resulting discrete theory is based on a consistent set of 
algebraic equations, but as one can see in equation 1)10(1 . one is enforcing the constraints with some variables evaluated 
at instant n and some at instant n + 1. For small stepsizes, this implies that the constraints with all the variables 
evaluated at the same instant of time are approximately preserved. The resulting theory therefore cannot be called 
mimetic in the traditional sense of the word, although it can do a good job of preserving (approximately) the discrete 
constraints. 

In the next two sections we apply this technique to linearized general relativity. We will see that the resulting 
theories do not determine the Lagrange multipliers, preserve the constraints automatically, and therefore are mimetic 
in the traditional sense of the word. We will not discuss the case of full general relativity here, but in several examples 
we have considered elsewhere for the non-linear theory (cosmologies , Gowdy spacetimes) the Lagrange multipliers 
are determined. Therefore it is unlikely that this method will yield a mimetic formulation for full GR. However, as 
we argued above, it will yield a formulation that approximates general relativity well in certain regimes and in such 
regimes the discrete constraints are enforced approximately very well. We believe it is likely that this is "as close as 
one will get" to a mimetic formulation of full general relativity. 



In this section we will apply the technique we described in the previous section to linearized general relativity 
written in terms of the traditional variables. We assume the background is the Minkowski metric. 



III. LINEARIZED GENERAL RELATIVITY IN TERMS OF METRIC VARIABLES 



A. 



Continuum formulation 



We start with the Arnowitt, Deser and Misner (ADM) form of the action of general relativity. 




(11) 
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where 



(12) 

C" = -2^"^^ (13) 
and the variables [qat, N, Na) are related to the four dimensional metric ^'^^g^u through, 

(14) 

-1/2 

.900 ) (15) 



qab = 


yabi 


N = 




Na = 





(16) 

The indices a,b,c run from 1 to 3. y/q is the determinant of the spatial metric qab and ^^^R is its Ricci curvature 
scalar. The momenta tt"'' are related to the extrinsic curvature of the space-like surfaces x*^ — t = constant through 
tt"'' = — [K"''' — q°'''K^] and indices are raised and lowered with the spatial metric. The semicolon denotes covariant 
differentiation with respect to the Christoffel connection of the spatial metric. Variation with respect to tt"'', qab, N, Na 
yields the Einstein equations. In particular variation of N, Na gives rise to four constraints C = 0, C" = usually 
referred to as (super)Hamiltonian and momentum (or diffeomorphism) constraints. 

We have chosen the ADM action since it is one of the most traditionally used in general relativity. Modern numerical 
implementations favor the use of formulations in which the evolution equations are manifestly symmetric-hyperbolic. 
This is not the case for the ADM equations. In principle there is no obstruction in applying our technique to any action, 
but it just is the case that there has been little investigations about formulating the symmetric-hyperbolic formulations 
as deriving from an action principle. This will require further study and therefore we decided to concentrate on this 
paper on the ADM action for simplicity. 

We now consider that the spacetime metric is given by a static background metric plus small perturbations g^i^ — 

^'^'^g^iJ + hfj_^. For simplicity we make the further choice that the foliation is such that the zeroth order shift Na = 
and the zeroth order extrinsic curvature is therefore zero tt"'' = 0. The constraint equations to leading order in the 
perturbations are given by Q , 



C 



ll-ab — n-.a 



hab^^'^R'^' 



(17) 
(18) 



In these expressions p°-^ is the linear portion of the canonical momentum tt"'' and R°-^ is the Ricci tensor of the 
background metric. The action for the linearized theory is. 



S 



p^^hab - N'^°^H - N^C - N'-^^C 



(19) 



where we have kept track of the order in the perturbation expansion of the lapse and the shift (recall that we assume 
zero shift in the background). The constraints C, C" are given by the expressions above, where only terms up to order 
linear have been kept. The quantity 



P^'pab-lp' 



- iKah^'' + 2h.,ah-''.,t + hh-',at - hhab^^^ R'"' 



(20) 



is a true Hamiltonian density (not a constraint) that is responsible for the evolution of the canonical variables, and is 
multiplied in the action times the lapse of the background space-time. 

At this point we can make an important observation. The momentum constraint (|17|l is only a function of the 
momenta p"'' (the covariant derivative is with respect to the background metric) and the Hamiltonian constraint 
(|18|l is only a function of the configuration variables g"''. Therefore, as we discussed in section II, our discretization 
technique will not determine the value of the Lagrange multipliers. The resulting theory therefore can only either be: 
mimetic or inconsistent. We will proceed to show that the resulting discrete theory is indeed consistent. 



5 



B. Discretization 



We start by discretizing the linearized action, S = ^^=1 -^("' ^ + where, 



L{n,n + 1) = XI 5Z {Pab{n,m){hab{n + l,m) - habin,m)) 

fa \a,b— 1 

~N{n, rh) [habin, rh + Ca + Sb) - hab{n, m - Ca + Bb) 

-hab{n, m + Ca- Cb) + hab{n, rh - Ca- Cb) 

- haa{n, m + 2eb) + 2haa{n, m) - haa{n, rh - 2eb)] 

-Na{n, m) [2pab{n, rh + Cb) - 2pab{n, rh - Cb)]} - H{n, rh) 



(21) 



and. 



H{n,rh) = ^ 



a,6=l 



Pab{n,rhf - -Paa{n,fh)pbb{n,fh) 



(22) 



1 



1 2 

- {hab{n, rh + ec)- hab{n, rh - Bc)) 



a.b.c—l 

- {hab{n, rh + Cc) - hab{n, rh. - Cc)) {hac{n, rh + et) - hac{n, rh. - Cb)) 
(haain, rh. + e,,) - haa{n, rh - Cc)) {hbtin. rh. + Cc) - hbb{n, rh - Cc)) 

+2 {hcc{n, rh + Ca) - hcc{n. rh - {hab{n, rh. + Cb) - hab{n, rh - eb)) 
+hcc{n, rh) [hab{n, rh + Ca + 4) - hab{n, rh + Ca - Cb) 

+hab{n, rh - Ca- eb) - hab{n, rh - Ca + eb) 

where we have assumed that the background metric is Minkowski and wc have chosen the zcroth order lapse equal 
to unity and we have dropped the (1) superscript from the first order lapse and shift. We have also chosen a 
centered prescription for spatial derivatives, with the following conventions, i.e. 4>{i)^x = + 1) — — 1) and 
4'{i),xx = (^{i + 2) + 4'{i — 2) — 2(/)(i) and similarly for higher derivatives. This choice of prescription is needed for two 
reasons: i) it ensures that "summation by parts" (ignoring boundaries) is satisfied, which is important when taking 
variations of the action; ii) it makes the successive application of two first derivatives the second derivative, etc. This 
is important when proving mimetism. 

The Lagrangian is the generator of the canonical transformation that materializes evolution from instant n to 
instant n + 1. Specifically, we will introduce the canonically conjugate momenta as we discussed in the previous 
section, 



-PaV"- + 1, ™) = Pab{n, fh) 

Pl^{n + l,rh)^Q 
P^(n + l,m) = 
P<f(n + l,m) = 

Pab{n, rh) = Pab{n, m) + N{n, m-Ca- Cb)- N{n, rh - ea + eb) 
-N{n, m + ea- Bb) + N{n, rh + ea + eb) 

3 

-Sab {N{n, m - 2ec) - 2N{n, rh) + N{n, m + 2ec)) 



(23) 
(24) 
(25) 
(26) 
(27) 



c=l 
3 



+ ^ V] [{hab{n, rh) - hab{n, fh - 2^)) - {ha,b{n, rh + 2ec) - hab{n, fh))] 
^ c=l 

1 ^ 

~ X] [(^ac(?^, fh + Bb-Bc)- hac{n, fh - eb- 4)) - {Kc{n, fh + eb + Cc)- hac{n, fh - eb + 4))] 



c=l 
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^ X! nebcin, rfi + ea- Be) - hbc{n, m - Ca - Cc)) - {hbc{n, m + Ca + ec) - hbc{n, m- ea + Cc))] 

c=l 

1 ^ 

^<5afc ^ m) - /id<i('^, m - 2ec)) - {hdd{n, m + 2ec) - m))] 



2 

c,ci=l 

+ -|- ^ [/icd(n-, m + - Cc) - /icd(n, m - Cd - Cc) - /icdCj^, + + e^) + hcd{n, m + e c - Cd)] 

c,d=l 

1 ^ 

+ 2 X! [('^cc(?^, m + 4 - eh) - hcc{n, m - Ca - Cb)) - {Kdn, rn + 4 + Cb) - hcc{n, m - Ca + Cb))] , 

P^j(n, m) = - [habin + 1, to) - hab{n, m)) + 2pab{n, rh) (28) 

3 

~ ^Pcc(n, rii)5ab + Na{n,rn- e^) - Na{n,rh + Cb) + Nb{n, rh - Ca) ~ Nb{n, rh + 4) 

c=l 

3 

P^" (n, to) = ^ [2pQ6(n, TO + efc) - 2pab{n, rh ~ Cb)] (29) 

b=l 
3 

P^{n,m,)= ^ [/iab(n, TO + 4 + efc) - ft,ab(n,TO - 4 + Cfc) (30) 

a,fc=l 

-hab{n, rh + ea- eb) + hab{n, rh - ea - Cb) - haa{n, rh, + 2e{,) + 2haa{n, rh) - haa{n, rh - 2e{,)] 

The system has four primary constraints H24I27|I . Preserving these constraints in time imphes, via H28I30|1 that the 
linearized Hamiltonian and momentum constraints are satisfied, 

3 

= 2 ^ [P':,{n, rh + e-fc) - P^U", ^ " ^b)] = (31) 

h=\ 
3 

C = ^ \hab[ri, TO + Ca + efc) - /la6('T-, TO - 6 q + Bf,) - ft.ofc(n, TO + Bq - Bfc) + ft.ofc(f^, TO - Ba - Bfc) 
a, 6=1 

-haa{n, rh + 2ef,) + 2haa{n, rh) - haa{n, rh ~ 2ef,)] = (32) 

Constraints (|24I27|) can be imposed strongly, the second constraint determines the variable Pab- This eliminates the 
variable pab and its canonically conjugate momenta from the theory. 
We now combine (|23|l and H27I) to get the evolution equation for P'\ 

P,''fc(n + l,TO) = PiUn,m) (33) 
—N{n, rh — Ca — eb) + N{n, rh — Ca + eb) + N{n, to + Bq — efc) — N{n, rh + Ca + Bfc) 

3 

+Sab ^ iN{n, rh - 24) - 2N{n, rh) + N{n, rh + 24)) 

c=l 
3 



^ [{hab{n, rh) - hab{n, rh - 24)) - {hab{n, rh + 24) - hab{n, rh))] 

c=l 
1 ^ 

+ - ^ [[hac{n, rh + Sb ~ Cc) - hac{n, rh ~ eb - 4)) - [hadn, to + 4 + 4) - /iac(", m - 4 + 4))] 

1 ^ 

+ - ^ [(/i6c(", TO + 4 - 4) - /i6c(»^, TO - 4 - 4)) - (hbc{n, TO + 4 + 4) - hbc{n, rh - Ca + 4))] 

c=l 
1 ^ 

+ -46 ^ [(/id<i(»^, m) - hddin, rh - 24)) - {hdd{n, rh + 24) - hdd{n, rh))] 



2 

c.d=l 



6 ^ 

^ [{hcd{n,rh + 4 - 4) - hcd(n,rh- 4 - 4)) - {hcd{n,rh + 4 + 4) - hcd{n,fh + 4 - 4))] 



2 

c,d=l 
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1 ^ 

^ 2 ^ [(^ccCj^,™ + Ca - ei,) - /icc(f^, m - Ca - eb)) - (/icc(n,m + Ca + eb) - /icc(", m - e a + Cb))] , 
and from (|28|l we get the evolution equation for h, 

3 

/iQf,(n + 1, m) = hab{n, m) + 2P^i,{n, m) - Sab Pffin, m) (34) 
+7VQ(ri,, m - ef,) - Na{n, m + eb) + Nb{n, m ~ Ca) - Nb{n, m + Ea) 

~2N{n, fn — Ca — Cb) + 2N(n, m — Ca + eb) + 2N(n, m + — ei,) — 2N{n, m + Ca + et) 

3 

- ^ [2hab{n, rn) - /ia6(n, m - 2ec) - /^a6(?^, m + 2ec)] 

3 

+ ^ [hac{n, m + eb- Cc) - hac{n, m - Cb - Cc) - hac{n, m + Cb + Cc) + hac{n, m - eb + ec)] 

c=l 
3 

+ ^ [^6c("-, m. + Ea - ec) - hbc{n, m ~ ea - ec) - hbdji, m + Gq + ec) + hbc[n, rn - ea + 4)] 

c=l 

1 ^ 

+ 2'^°'' X! [2hdd{n, rh) - /i<i<i(»^, ™ - 2ec) - hddin, rn + 2ec)] 

1 ^ 

"2'^"'' ^ [^cd{n, rn + ed- Cc) ~ hcd{n, m - e^ - Be) - hcd{n, rh + Cc + ed) + hcd{n, rn + e^ ~ e^)] 

3 

- X! ['^ccC"-, TO + ea - 4) - /icc(", rn - ea - eb) - hcc{n, rh + ea + Cb) + hcdn, rh - Sa + eb)] , 



A first point to be noted is that the evolution equations have resulted in an explicit evolution scheme. This is 
usually not the case, it is a particularity of the linearized theory that the evolution is explicit. It should be noted 
that the evolution equations obtained are just a straightforward discretization of the evolution equations one would 
obtain in the continuum by working out the variations of the continuum action. 

We have checked, using a computer algebra code, that the evolution equations H34I33() exactly preserve the con- 
straints H31I32(I . or more precisely that, 

Cain + l,m) = Ca{n,m), (35) 

3 

C(n-t-l,m) = C{n,m) +'Y[Ca{n,m + ea) - Ca{n,m- ea)]. (36) 

a=l 

This result was expected since we used differentiation operators that ensure that mixed discrete spatial derivatives 
commute, and that one can integrate by parts (more precisely "sum by parts"), and that is all that is needed in a 
linear theory on a Minkowski background to show that the constraints are preserved upon evolution. It is interesting 
to compare this result with that of Meier (Sj] . He finds a mimetic discretization of linearized general relativity around 
Minkowski spacetime, but using staggered grids. This is a natural approach, for instance, in electromagnetism and 
Yang-Mills theory (and it is the one we will take in the next section where we deal with gravity with the Ashtekar 
variables) . 

It would be interesting to generalize these results to the case of linearization around a static background. In that 
case it is not obvious that the formulation would result automatically mimetic. In fact, the failure of the Leibnitz 
rule at a discrete level implies that it will be difficult to find a mimetic formulation since the equations now will have 
non-constant coefficients and one will need Leibnitz' rule to show conservation. Our formalism will yield a consistent 
formulation, but it is possible that it will require determining the Lagrange multipliers. 

C. Stability 



We have discretized the time derivatives without centering them (that is, we have used a stencil that is first order 
accurate only). The reason for this is that the canonical theory is much cleaner with only two levels in time involved 
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in the derivatives. It is possible to use derivatives that are second order accurate in time and use our construction 
by rewriting the theory in terms new variables in such a way that the resulting theory has derivatives that are first 
order accurate, but we will not do this here. 

The spatial derivatives, on the other hand, were centered (this was required in order to have summation by parts). 
The resulting scheme is therefore "forward in time centered in space" , a recipe that is not stable, for instance, for 
the advection or the wave equation. We therefore would like to check if our scheme is stable. To simplify things, we 
will consider H34I33|I and make the following assumptions: the metric and extrinsic curvatures are diagonal and only 
depend on the coordinates t, x, the lapse is unity and the shift is zero. The resulting equations therefore are, 

1 ^ 

Piiin + l,m) = P^^{n,m) - [2ft,cc(n, m) - hcc{n,Tn- 2ei) - hccin,m + 2ei)] , (37) 

P^'^in + l,m)^ PiUn, ™) + ^ [2/^33(^^, m) - /i33(n, m - 2ei) - hsain, rh + 2ei)] (38) 

Pssi^ + 1, ?n) = Pssin, to) + ^ [2/122 (n, m) - /i22(n-, rh - 2ei) - ft.22(?^, to + 2ei)] (39) 



hii{n + l,m) = hii{n,m) + 2P^^{n,m) - ^P^f(n,m) (40) 

/=i 

1 ^ 

^ [2hdd{n,m) - hdd{n,m- 2ei) - hdd{n,m + 2ei)] 

d=2 



h22{n + l,m) = h22{n,m) + 2P^2in^m) -^P^j{n,m) (41) 

/=i 

- 2 [2/^22 ('^i to) - h22{n, TO - 2ei) - /i22(", m + 2ei)] 
+ i [2/i33(n, rh) - /i33(", rh - 2ei) - /i33(", rh + 2e2)] 

3 

h33{n + l,m) = h33{n,m) +2P^s{n,m) -J^Pffi^^m) (42) 

/=i 

-- [2/i33(n, rh) - h^zin, rh - 2ei) - h33{n, rh + 2ei)] 

+ i [2h22{n, rh) ~ h22{n, rh - 2ei) - /i22(n, to + 2e2)] 

As a test case, we concentrate on a subfamily of solutions of the equations, in which hn — Pn — and ft-22 — — /133 
and P22 = —^^33. In that case, the equations reduce to, 

P2''2('T' + 1,to) = P22{n,rh) - i [2/i22(n, to) - ft,22(«, to - 2ei) - /i22(n, to + 2ei)] , (43) 

/i22(n + 1,to) = h22{n, rh) + 2P22in,rh) - [2h22in,rh) ~ ft,22(n,TO - 2ei) - ft,22(", to + 2ei)] . (44) 

We have performed a Von Neumann analysis of this system and confirmed that the scheme is stable provided the 
Courant factor is less than one. So at least for this particular subcase the scheme is stable. A more complete analysis 
is needed to guarantee stability in general. 
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IV. LINEARIZED GENERAL RELATIVITY IN TERMS OF ASHTEKAR VARIABLES 

A. Continuum formulation 

We will now apply the technique we outlined in the previous section to general relativity linearized around Minkowski 
space using the Ashtekar formulation. The formulation of linearized gravity with the new variables was first discussed 
by Ashtekar and Lee Q • The discussion presented in that paper required the use of complex variables if one was to 
describe general relativity with metrics with a Lorentzian signature (alternatively, one could consider real variables, 
but then the theory described the Euclidean signature sector.) Developments that have taken place in the field since 
the publication of that paper that allow to consider the Lorentzian sector using real variables |lCt lll| , but we will see 
that the discretized theory is more problematic in this case and we will not discuss it in detail in this paper. 

The Ashtekar canonical variables consist of a set of triads with density weight 1, and a (complex) 5*0(3) 
connection Aai- In this notation a,b, . . . are spatial vector indices and . . . range from 1 to 3. Following Ashtekar 
and Lee we omit using tildes to denote density weights since in this context they do not play an important role. To 
linearize the theory around Minkowski, we choose a fixed background (£^°' — Eg'^jAai — 0) in the phase space and 
consider fluctuations around it. In Cartesian coordinates, E'q' = (5°*. The triad is therefore given by, 

£;"' = (5" + e", (45) 

and therefore the background metric has components q"''' — (5°^ and its determinant is unity and therefore density 
weights are all trivial. We will denote by Aai the fluctuations of the connection. The Poisson bracket of the canonical 
variables is {e'^^{x), Abj{y)} — iS^6jS{x — y). 

The Ashtekar formulation has, in addition to the usual diffeomorphism and Hamiltonian constraints of the metric 
canonical formulation of general relativity, a set of additional constraints that make the formulation invariant under 
triad rotations. The additional constraints take the form of a Gauss law, which linearized will read, 

Gl = dac'^' + e'^'^Aaj = 0, (46) 

where from now on the subscript L means we have kept the minimum required number of terms in the perturbative 
expansion. In spite of the second term, one can check that if one computes the Poisson bracket of two Gauss laws, they 
commute, that is, they form an Abelian algebra. The internal symmetry group of the linearized theory is therefore 

Ignoring boundary terms, the (super)IIamiltonian for general relativity can be written as. 



H = I (fx 



(47) 



and in the full theory it vanishes identically The parameter (3 is called the Immirzi parameter and the parameter a 
is equal to +1 for the Euclidean case and —1 for the Lorentzian signature. Classically, different values of the Immirzi 
parameter correspond to different representations of the same theory. The quantities are the spin connections 
compatible with the triads, defined by, 

d[aEl^+e]kTm^O, (48) 

where the iJ's are the triads (without density weight), related to the Ashtekar variables by Ef = det{E)Ef , or 
equivalently, Ef — E^E^e'^^'^eijk- Indices are lowered and raised with the flat Euclidean metric. One can obtain an 
explicit expression for the spin connection in terms of the triads, 

F^ = ef(9,4-9„Tr(e)<5^). (49) 

To study it in the linearized theory, we need to choose a lapse and a shift. The natural choice is to use as zeroth 
order lapse and shifts the ones that would preserve the spatial background metric explicitly time- independent. This 
corresponds to a lapse iV = 1 and a shift N°' = 0. So we will write Nl = \ -\- v and iV£ = and these will become 
Lagrange multipliers in the linearized theory. The super-Hamiltonian then separates into two pieces, one that acts as 
a Hamiltonian and another piece that is given by the Lagrange multipliers times constraints of the linearized theory. 
These constraints are, 

C^. = -»/a6 = 0, (50) 

= -*ef/a^fc = 0, (51) 
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where /^^ = 29[aA^j is the hnearized field strength, and the first one is the linearized momentum constraint and the 
latter the linearized Hamiltonian constraint. The non- vanishing Hamiltonian for the linearized theory is given by, 



= 1 d^x {2etf^,et + {A:a1, - A'^A^) (52) 

[(r^ - <Ji3A:){rl - a/34)) - (r^ - '^/34)(r? - ^/3A?)] 



f3' 



/32 



B. Discretizing the full theory on the lattice 



In this section we review some results of reference |^|5j where we discretized general relativity on the lattice. In the 
next section we will particularize these results to the case of linearized general relativity. We start by considering an 
action for general relativity written in terms of Ashtekar's variables (see for instance 12] and the book by Ashtekar 
M page 47), 

''e'^^FI^^H (53) 

where N and are the lapse and shift and H the supcr-Hamiltonian H47|) . We will particularize to the Euclidean 
case a — 1 and choose the Immirzi parameter 13—1 which correspond to the original form of Ashtekar's variables, 
for simplicity (see section VI for more details). From now on we will not assume Einstein's summation convention 
and present the summations explicitly, since many expressions would otherwise be confusing. The Lagrangian can be 
discretized as follows, 

^n,v\"'n,v ""n,v) / , ^^n.v \"'n.v 



^ n^v\ n,v n^vJ /_ ^ n.v\ n.v n,D/ 

a a,b 



L{n,n+l) = 

V 

+ ^aa,„,. - l) +f3n,. {hl^ {hl^)^ - l) 



(54) 



where /i"^]^ represents an holonomy along the a direction at instant n + 1, h^^ represents the "vertical" (time-like) 
holonomy. The holonomy associated with a plaquette in the a(3 {a ^ (3) plane (a, /3 = . . . 3) is 



and 



We will assume that the holonomies are matrices of the form h — h^T^ where = I and T"^ — —ia"" where ct" 
are the Pauli matrices. The indices n, v represent a label for "time" n and a spatial label for the vertices of the lattice 
V. The elementary unit vectors along the spatial directions are labeled as ea, so for instance n + ei labels the nearest 
neighbor to n along the ei direction. The unit vector in the timelike direction is eo and we chose /i" ^^eo = ^n+i v 
The quantities ^ are elements of the algebra of su(2) and a and /3 are Lagrange multipliers, the last two terms 
of the Lagrangian enforcing the condition that the holonomies are elements of SU{2). We use the usual conventions 
of lattice gauge theories in which one has oriented links and the natural variables are the holonomies in a given 
orientation and based at a given vertex. If we need to traverse back, as in the case of closed loops one then considers 
the adjoint of the holonomy based at the vertex one is ending at. 
The discretization of the field tensor is based on, 

Kb^-\Tr[Kv-hn,.)r] ■ (57) 

Instead of working out the equations of motion for this action, we will, in the next section, particularize it to the 
linearized case and work out the relevant equations of motion, which is equivalent to working the equations first and 
then linearizing if appropriate perturbative orders are kept. 
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C. The linearized theory on the lattice 



We now proceed to linearize the action. We start with the holonomies. The exphcit form of the hnearized holonomy 



(58) 



where we have dropped the subscript n we used in the last section to indicate the time level in order to make the 
notation more compact (but we will make it explicit when things are evaluated at n + 1). In this equation (/)"*T* 
is an element of the algebra that can be viewed as a "phase" (it corresponds to the logarithm of the path-ordered 
exponential of the connection along the direction a). The holonomy of a plaquette in the plane a/3 is (neglecting 
higher order terms), 



i i 

The first order contribution is, 



^afJk ^ 
Iv — 



A0k 



and the second order contribution is. 



a/3 ij _ 



2v 
2v 



(59) 

(60) 

(61) 
(62) 



We now linearize the expression for K defined in the previous section, by noting that to first order. 



(63) 



and ignoring higher order terms we get that 



(64) 



with 



T^abk — ^afcA:,, , "^/,,ar6fc ^,&rafc\ , \ ^ ( ^ai ihk hi iak\ 



(65) 



We now consider the first term in the discretized Lagrangian (|54|l . Substituting the expression for the holonomy 
around a plaquette H59|l we get the following identity, valid up to second order. 



iTr 
4 



Epa //jflO _ 1,0a \ 



(66) 



ak 



and we note that when one considers the sum over all vertices, the first term of the right hand side yields a total 
derivative with respect to time that can be ignored in the Lagrangian. 

For the second term in (|54|) we use (|59ll and (|64|l . getting the following identity, valid up to second order. 



^ rp r \ ^ T^ab ( udb uba \l \ ^ f abk ^abk i ^abk^abk i T^abk^ahk\ 



(67) 



ab 



abk 



When one considers the sum over all vertices the first term on the right hand side of this expression vanishes. The 
resulting Lagrangian therefore can be written as. 



L = - 



EE' 



■V \'rn+l,v 



'^<^tjk{<l)'J(f>'n+l^v 



ijk 
2v 



ijk 



aijk 



ab 



ij k ij k 



(68) 
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Now that we have an expHcit expression for the Lagrangian we can proceed to identify the various terms. The 
theory has the following Lagrange multipliers: 0°' the "vertical component of the phase" (which plays a role analogous 
to the time component of the vector potential in Maxwell theory) and the linearized lapse and shift. These quantities 
multiply times the constraints of the linearized theory. Explicitly, the momentum and Hamiltonian constraint read, 

b 

a = 5:e.,fc$^^. (70) 

ijk 

In order to get Gauss' law, we first take the variation of the Lagrangian with respect to the Lagrange multiplier 
to get, 

C^l = ^( + ^n+l,v ) = (71) 

a 

where 

k 

C+1,. ^ -<5" - C:.-e. + E (C-e. + C + C-e. + C+l,«-e„ ) (73) 

k 

At the moment this does not appear to be a true constraint since it involves variables at instant n and at instant 
n + 1. To see that it actually is a constraint, we will call ^n+i v the component in the direction ia of a quantity that 
we will think of as an "electric field" (in the sense that it is the quantity that satisfies the usual form of Gauss law) 
and we will call „ the component in the direction —ia, both at point {n + 1,1;). To make this more transparent, 
we need to see how they transform under gauge transformations. To leading order the field ej^ „ is eg „ = J^i 6°'^T'^ ■ 
We then define 



n,v ~ ^n,v ' A"-n,v '^0,n,v \"'n,v ) "■n,v '^Q,n,v \""a,v ) J ~ ^n.v ' / , "^ajfc^lu J- i \'^) 



4' 

jk 



with the second equality valid up to second order. By inspection one sees that the field e ° ^ is an element of the 
algebra that transforms like an electric field at (n, v) under gauge transformation. One can also show the following 
identities, valid to first order, 

Cn+l,t) = ^ ^ ^n+l.v ^ ~ i^n.v)^ ^n,v ^n,v (75) 
i 

in+l,v = Cra+1,D ^ — ~(^n+l,t)-e„)^ Cn+l,t)-e„ ^n+l,D-ea (76) 

i 

from which one immediately sees that the quantities we identified as components of the electric field have the appro- 
priate transformation properties under gauge transformations. 

Therefore we can identify H71|) as the usual intuitive expression of Gauss' law stating that field lines cannot emanate 
from a point in vacuum. 

We now turn our attention to the equations of motion. Given the Lagrangian (|68|l we work out the equations of 
motion from the canonical transformation. We start by computing the canonical conjugate momentum to e, 

p{e)ak _ dL{n, n+1) _ 

pi-yi^ _ _ dL{n^n+l) ^ _^,ofc ^ 2 ^ e^'^^ cS>^ = 0. (78) 

Therefore the dynamics of P^'^-' is trivial. However, the last equation can be viewed as an evolution equation for 
through (EOJ), 

rn\i,. = - <t^f + <t^f+e. - 2 E • (79) 
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Notice that by adding over indices a belonging to a given plaquette equation H79() . one effectively gets an evolution 
equation for all the horizontal 0's in the plaquette. This is due to the fact that the vertical contributions in (|tiO|) will 
cancel out in pairs when adding through the plaquette. Explicitly, 

where in the last term Vd is the vertex in which the link d originates and Pab is the plaquette spanned by a and b. 

We now consider the momentum canonically conjugate to (j). We start by computing the canonical conjugate 
momentum at instant n + 1 

p(</,)afc _ dL{n,n+ 1) ^ak , ST' akz / ,ai , ,0i , ,0i N (oT\ 

The momentum can be written in terms of the electric field in an expression that sees parallels the usual relation 
between the electric field and the canonical momentum in the lattice ^""^ as, 

^i+M = ^ak - inXl,v + E ^"'^^ ^n+l,v (82) 

i 

and in terms of it, Gauss' law l|71|l can be written as, 

f-^k _ _ ( r>Wak _ r,{(i))ak _ aki ( lai 

a \ i 



(83) 



This final expression for Gauss' law is a genuine constraint, in the sense that all variables are expressed at the same 
instant of time. 

We now compute the momentum conjugate to (j) ^t instant n. 



p(4,)ab _ dL{n,n+l) _ ^ + S .S" ( -u' + ly' -] 

i 

+ E ^-br {-2V, + 2i^„_e, + 4>v - €+e^ - + '^^-e, ~ '^^ + Se! 

i 

-2e^ + 2et,^ + C+i,.) + 2<5ab ^ '-^i " 



(84) 



One still needs to replace the expressions for the e's and for the 0's evaluated at instant n + 1. The resulting 
substitutions lead to lengthy expressions that are not particularly illuminating, and will not be needed in what 
follows, so we will not display them here. We point out however, that the resulting scheme is not an explicit one 
for the P^-^^'s. Since these variables do not arise in the constraints, we do not need this evolution equation to show 
mimetism. 

One now needs to show that the evolution is mimetic, that is, it preserves the discrete constraints (|69l 1701 183II . 
Using the evolution equation 18U|) one gets that, 

^n+l,v ~ + ~ Ctj+e^ , (85) 

- a + 4 E(C+e„ - O- (86) 

a 

To study the time evolution of Gauss' law one needs equations (|79I81I84|I and one gets that 

Gn+1,D — ^ri,!) ■ (87) 

As in the previous section, we have checked these identities using computer algebra. 
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V. DISCUSSION AND CONCLUSIONS 

The consistent discretization scheme is such that it yields a set of discrete equations for the evolution equations 
and the constraints of general relativity that is compatible, that is. all can be solved simultaneously It docs so at the 
price of determining the Lagrange multipliers (the lapse and the shift). In the linearized case we have shown that one 
can discretize the theory in such a way that the Lagrange multipliers are not determined and nevertheless the theory 
is consistent. 

When one discretizes a theory there is always an ambiguity in how to proceed. Among the ambiguities we have the 
dependence on how one chooses to represent the derivatives. What we have found is that in the linearized case one 
can choose certain derivative operators for which the Lagrange multipliers are not determined. It should be noted 
that the consistent discretization scheme would work even if one did not choose the derivatives this way, but the 
Lagrange multipliers will be determined in order to have a consistent set of equations. This is true both in the case 
of the traditional variables and also in the Ashtekar variables. In the latter case there is an additional element that 
is the presence of an extra constraint: the Gauss law. We have also chosen a specific way of discretizing the theory 
in such a way that Gauss' law is implemented exactly in the discrete theory (this is standard in Yang Mills theory 
on the lattice, and implies that the discrete formulation is gauge invariant, and also that these discretizations are 
mimetic, though this is rarely emphasized in the Yang-Mills literature). In the case of the traditional variables, one 
can also associate mimetism with gauge invariance. The action we chose to work with is invariant under linearized 
coordinate transformations of the form h'^^^ — h^i, + S,{^^v) ■ The discrete action, if one chooses a derivative operator 
such that the second derivatives coincide with the derivative of a first derivative and satisfies summation by parts, 
is invariant under a discrete version of the above symmetry. This symmetry is generated canonically by the discrete 
constraints. This explains in a geometrically nice way why mimetism was possible in the linearized case. 

In the case of Lorcntzian general relativity written in terms of Ashtekar's variables, the presence of the terms 
(r^ — (j(3A^g) in the Hamiltonian make it more difficult to discretize the action in such a way that the Gauss law is 
preserved exactly. This is because the F^'s have to be written in terms of the triads and the resulting expressions are 
not easy to discretize on the lattice preserving the internal symmetry (unlike the A\^s which are readily discretized 
by considering a parallel transport operator along the elementary links). It may be possible using dual lattices to 
discretize such terms in an invariant way, but this will require further study. There is no problem applying our 
discretization technique to this case directly, but what will happen is that the internal symmetry will be broken, the 
Lagrange multipliers associated with the Gauss law will be determined and the resulting discretization will not be 
mimetic in the traditional sense of the word. It is clear that further work is needed before this kind of discretizations 
will be useful numerically. 

In the full nonlinear case either with the traditional or the new variables, the constraints involve both coordinates 
and momenta and therefore the application of our technique will determine the Lagrange multipliers and will therefore 
not furnish a mimetic discretization in the traditional sense. The resulting discrete theory is consistent, but it does 
so at the expense of determining the Lagrange multipliers. Based on what we learned from the linearized case, we 
can conclude that the only remaining possibility for a formulation that is mimetic in the traditional sense would be 
to implement the symmetries implied by the constraints exactly in the lattice. Since the symmetries implied by the 
diffeomorphism constraint are broken by the introduction of the lattice it appears unlikely that such a formulation 
would ever be found. 

The conclusion we can draw from this is that for the case of full nonlinear general relativity, the closest one can come 
to a formulation that preserves the constraints under evolution is the proposal of consistent discretizations we have 
introduced. Such proposal is not mimetic in the traditional sense in that it imposes the constraints by determining the 
Lagrange multipliers. This proposal has many new aspects that are currently in investigation. It has been successfully 
applied in cosmological examples and is now being studied in detail for the Gowdy space-times. If it works for this 
example, it is likely that it could be applied successfully in general, but this obviously requires further study. 

Something to be noticed is that it is not clear that the formulations we presented are going to be useful numerically. 
In particular, the fact that they are not based on manifestly hyperbolic equations. We have presented a first step 
towards showing stability of the scheme in a particular situation, but a fuller analysis should be carried out to 
determine if the scheme is stable in general. Numerical relativity codes also use more sophisticated time stepping 
techniques than the one we use. It is clear, however, that our method can accommodate more elaborate discretizations 
of the time derivatives and the calculations in this paper could be repeated in that case. Another interesting point 
would be to attempt to apply the techniques in this paper to the several manifestly hyperbolic formulations of the 
Einstein equations that have been proposed in the last few years. Unfortunately, few of them have been worked out 
in the context of an action principle, but this difficulty could presumably be remedied. This would also allow to study 
within our framework manifolds with boundaries. 

Another element of interest is the impact of the choice of the derivative operators on the construction of consistent 
discrete theories. The consistent technique will work no matter what derivative operators are chosen. But here we 
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have learnt that one can choose them in such a way that the Hnearized theory is automatically mimetic. We would 
like to argue that the level of accuracy with which the consistent discretization enforces the constraints is improved 
when one chooses a formulation that is mimetic at the linear level, at least for weak fields. The argument is simple. 
In the consistent discretization scheme the constraints that are enforced exactly have the form (j){qn,Pn+i) — 0. The 
constraints one would like to see enforced arc of the form (f){qmPn) = 0. Starting from the former, and using the 
equations of motion one has that (/)(g„,p„ + 0{p^)) = where the terms that correct pn are of order p^ (or 
or mixed but quadratic). This is true if the theory is mimetic in the linearized level. Otherwise one would have 
4>{ln,Pn + 0{p)) = 0. Therefore choosing a discretization that is mimetic at the linearized level, at least for weak 
fields, implies that the constraints are tracked more accurately in the full non-linear theory when one discretizes 
consistently. 

Summarizing, we have shown that the consistent discretization scheme we have introduced recently, when applied to 
general relativity discretized around Minkowski spacetime yields a formulation that is mimetic. That is, a formulation 
in which the discrete constraints are exactly preserved upon discrete evolution. We have also argued that for the full 
nonlinear case, the use of consistent discretizations appears as a possibility to yield a formulation that is close to the 
intention of mimetic formulations, although only approximately. 



VI. ACKNOWLEDGMENTS 



We wish to thank Manuel Tiglio for discussions and Luis Lehner and Olivier Sarbach for comments on the 
manuscript. This work was supported by grant nsf-phy0244335, NASA-NAG5-13430, DID-USB grant GID-30, Fonacit 
grant G-2001000712 and funds from the Horace Hearne Jr. Laboratory for Theoretical Physics and the Abdus Salam 
International Center for Theoretical Physics. 



A good guide to the mimetic literature can be found in |http://www.matli.unm.edu/'stanly/mimetic.html| 

M. Miller (private communication); Also see the talk by J. Frauendiener, 

http : // online. kitp .ucsb . edu/ online / gravity 03 /frauendiener / 

D. L. Meier, Astroph. J., 595, 980 (2003) arXiv:astro-ph/0312052 . 

C. Di Bartolo, R. Gambini and J. PuUin, Class. Quant. Grav. 19, 5275 (2002) |arXiv :gr-qc/0205123| . 
R. Gambini and J. PuUin, Phys. Rev. Lett. 90, 021301 (2003) arXiv:gr-qc/0206055\ 
R. Gambini and J. PuUin, Class. Quant. Grav. 20, 3341 (2003) arXiv:gr-qc/0212033 . 

R. Arnowitt, S. Deser, C. Misner, in "Gravitation: An introduction to current research", L. Witten, editor, Wiley, New 
York (1962). 

V. Moncrief, Ann. Phys. (NY) 88, 323 (1974). 

A. Ashtekar and J. H. Lee, Int. J. Mod. Phys. D 3, 675 (1994). 

J. F. Barbero Phys. Rev. D 51, 5507 (1995) arXiv:gr-qc/94100T4l. 

T. Thiemann, Class. Quant. Grav. 15, 839 (1998) arXiv:gr-qc/9606089 . 

R. Capovilla, T. Jacobson, J. Dell and L. Mason, Class. Quant. Grav. 8 (1991) 41. 

A. Ashtekar, "Lectures on non-perturbative quantum gravity" World Scientific, Singapore (1991). 



